Equilibrium long-ranged charge correlations at the surface 
of a conductor coupled to the electromagnetic radiation II 
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Results of a previous article with the same title are retrieved by a different method. A one- 
component plasma is bounded by a plane surface. The plasma is fully coupled to the electromagnetic 
field, therefore the charge correlations are retarded. The quantum correlation function of the surface 
charge densities, at times different by t, at asymptotical large distances R, at inverse temperature j3, 
decays as — l/(87r 2 /37? 3 ), a surprisingly simple result: the decay is independent of Planck's constant 
h and of the time difference t. The present paper is based on the analysis of the collective vibration 
modes of the system. 
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I. INTRODUCTION 

In a previous paper we studied the asymptotic 
form of the two-point correlation function of the surface 
charge densities on a plane wall bounding a conductor [in 
the special case when the conductor is a one-component 
plasma (OCP), also called jellium], taking into account 
the retardation and the quantum nature of both the one- 
component plasma and the radiation. The novelty was 
the retardation; instead of assuming as an interaction the 
Coulomb potential only, in [l[ the coupling was through 
the full electromagnetic radiation. In other words instead 
of assuming the velocity of light c to be infinite, the full 
Maxwell equations were used. 

The previous paper [l[ used the elaborate formalism 
of Rytov @, presented also in This formalism is 
macroscopic, using frequency-dependent dielectric func- 
tions. In the present paper, we retrieve the same results 
by a simpler method, based on the analysis of the col- 
lective vibration modes of the system. This method is 
partially microscopic. It has already been used in the 
non-retarded case [J], Sec. 4. 

We use Gaussian units. The OCP is made of point- 
particles of charge e, mass m, and number density n, im- 
mersed in a uniform neutralizing background of charge 
density — ne. We recall the geometry (Fig. 1). We use 
Cartesian coordinates; a point is r = (x, y, z). The OCP 
occupies the half-space Ai = {x > 0}, the half-space 
A2 = {x < 0} is vacuum; the two half-spaces are sepa- 
rated by an plane wall, impenetrable to the jellium, at 
x = 0. A point on the wall is R = (y, z). 

After long calculations, a very simple result was found 
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where j3 is the inverse temperature, (• • • ) T represents a 
truncated statistical average, a(t, R) is the surface charge 
density at time t and at point R on the surface. This 
value ([1]) surprisingly is independent of t and ft, boiling 
down to the classical result at time difference zero with- 
out retardation. 

We define the Fourier transform of a function /(R) as 



/(q) = J rf 2 i?exp(zq-R)/(R). 



(2) 



A result equivalent to |T]) is that the Fourier transform of 
its lhs, j3S(t, q), has a kink singularity at q = 0, behaving 
at small q like q/(4w). 

The present paper is organized as follows. Sec. II is a 
general exposition of the formalism of vibration modes. 
Sec. Ill describes the contribution of surface modes 
(which are localized on both sides of the wall). Sec. IV 
describes the contribution of transverse modes delocal- 
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FIG. 1: The geometry. 



2 



ized on the vacuum side. Sec. V describes the contri- 
bution of transverse modes delocalized on both sides. In 
Sec. VI, we recall the contribution of the longitudinal 
modes. Sec. VII is the Conclusion. 



II. COLLECTIVE VIBRATION MODES 

We consider a fluctuation of the surface charge density 
(7, of wave number q and frequency lo, such that a is of 
the form 



where the frequency-dependent dielectric function is 



er(R, t) — a^{t) exp(iq • R) + c.c. 



(3) 



where c.c. means complex conjugate. cr qLJ (i) is a com- 
plex quantity, vibrating at frequency to. The emitted 
radiation, of frequency u>, is described by the Maxwell 
equations (we use the microscopic ones, involving only 
the electric field E and the magnetic field B, averaged 
in a suitable way [5J). The charge density is p, the elec- 
tric current density is J, the velocity of light is c. These 
Maxwell equations are 



1 <9E 4tt , 

c at c 



V x E = — 

V • E 
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~c~dt' 
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In the quantum case, the quantities appearing in these 
equations are operators. In region A2, p = and J = 0, 
of course. 

There are solutions to these Maxwell equations which 
are superpositions of transverse waves, which we shall 
study first. In region Ai, for these transverse waves, p — 
0. Furthermore, the Maxwell equations (|4][7]) have to be 
supplemented by 



E 



4tt dJ 



(8) 



where uj p = (Anne 2 /m) 1 / 2 is the plasma frequency. This 
equation J5J) is obtained from assuming that the OCP can 
be treated in a hydrodynamical linearized approximation 
where the velocity is v(r, t) obeying the Newton equation 
mdv/dt = eE (the term ev x B in the Lorentz force is 
suppressed by the linearization) and the current density 
is J = env (the density is the constant n, again because 
of the linearization) . It should be remarked that in the 
equation l[8]). there is no damping term; this absence of 
damping term is valid for the small wavenumbers which 
will be considered, and is a property special to the OCP. 

Macroscopic Maxwell equations equivalent to (0][5|) can 
also be obtained. Combining ([4} and (8J, at frequency 
<x>, we obtain 
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(9) 



e(w) = 1 - 
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which is of the Drude form, without dissipation. Multi- 
plying (J6j> , where p = 0, by e does not change anything if 
e ^ (the case e — will be studied later). To an excel- 
lent approximation, the magnetic permeability is fi = 1. 
Thus we have obtained the macroscopic Maxwell equa- 
tions [f|, with now p = and J = 0. 

For these transverse waves, the wavenumber vector has 
components (fci, q) in region Ai and (k 2 ,q) in region A 2 , 
where k\ and k 2 are the x components (the wavenumbers 
have the same components parallel to the surface q — 
(QvtQz) as a consequence of the boundary conditions, as 
will be shown later). 

From the Maxwell equations (jU[7]) and (J8]) , the disper- 
sion relations are 



in region A2, and 



lo 2 = c 2 {q 2 + k 2 2 ) 



+ c 2 (g 2 + fc 2 ) 



(11) 



(12) 



in region Ai. Depending on the value of co, k\ or k 2 
are real or pure imaginary, as shown from (jlip and (|12p . 
If to 2 < (eg) 2 , both k\ and k 2 are pure imaginary. If 
(cq) 2 < uj 2 < to 2 + (eg) 2 , k\ is pure imaginary and k 2 is 
real. If oj 2 > uj 2 + (eg) 2 , both k\ and k 2 are real. 

In addition to these plane transverse waves, in region 
Ai, when the frequency is oj p (then e = 0) the Maxwell 
equations have solutions which are longitudinal waves 
Q . They have already been studied in i[ . 

is a 



The energy density in region A2 
1 



U 2 



-(E z + B 2 ); 



the energy density in region Ax is 
1 



8tt 
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A derivation of (|14[) is presented in the Appendix, as well 
as another form of it [6( . 

At x = 0, the existence of a surface charge density a 
is associated with a discontinuity of the x component of 
the electrical field 



Ana 



EZ 



(15) 



where the upperscripts + or — mean approaching the 
plane x — from the regions Ai or A2, respectively. The 
other conditions at x = are the continuity of the com- 
ponents of the electric and magnetic fields parallel to the 
surface. From the conditions at x = Q, one easily deduces 
that the (y, z) components of the wavenumbers indeed 
are q. Since the system is invariant by rotations around 
the x axis, general results can be obtained by choosing q 



along the y axis: q = (g y ,0); then g 2 
only modes which contribute to a. 



q . We consider 
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III. to 2 < (cq) 2 (POLARITONS) 

Since k\ and k 2 are pure imaginary, the solutions of the 
Maxwell equations are of the form (the subscript j = 1, 2 
denotes the region A x or A 2 ) Q 



Ej = [aj(t),bj(t),0] exp(iq y y - Kj\x\) 



Bj = [0, 0, dj(t)] exp(iq y y - Kj\x\) 



(16) 
(17) 



where, from (jlip and (TH?)) . cn 2 = \J (cq) 2 — lo 2 and 

cki = \Ju 2 + (cq) 2 — lu 2 . Since all fields are localized 

near the surface, these modes are called surface plasmons; 
polariton is also used when retardation, as here, is taken 
into account. 

From the Maxwell equations and the conditions that 
E y and B z are continuous at x = 0, in this Section there 
is a relation between lu and q: 



lu 2 = lu 2 /2 + (cq) 2 ^(lu 2 /2) 2 + (cqf 



(18) 



Therefore, in this Section, er qw will be called £r q . lu has 
its non-retarded value cu p /V2 only at large q. For the 
small values of q, that we are interested in, lu behaves as 
cq. From (fT8]) follows kik 2 = q 2 , a useful relation in the 
calculations (the detail of which is omitted) which follow. 

Taking into account (fT5")) and the continuity of E y and 
B z , the Maxwell equations ©, © [p = in ©] give the 
prefactors a,i,bi, di as functions of <7 q . One finds 



ai(t) = 

Kl + K 2 

= fe(t) = -»- 



-47rcr q (i), a 2 (t) = -- 



Kl + K 2 



d 1 (t)=d 2 (t) 



- k 2 
1 



-47rcr q (t), 

■47T(7 q (t). 



-4w(7 q (t), 

(19) 



CK2 Kl + K 2 

In region Ai, J can be obtained from ([4|). One finds 
Jia = - exp(ig a y - K X a;)a- q (t) + c.c, 
Ji y — — - exp(iq y y — Kix)&„(t) + c.c. (20) 

K 2 

From (UHJ) and (j2"0)) . one computes the energy densities 
(fT3|) and (fT4)) and the total energy iJ q . Later, we shall 
consider that the large area A of the wall goes to infin- 
ity; therefore, the oscillatory terms exp(±2iq y y) do not 
contribute to the integral on R. One finds 
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= AC^ |a q (i)| 2 + — |<7 q (t)| 



(21) 



One sees that i? q is the energy of a two-dimensional har- 
monic oscillator (two-dimensional because cr q (i) is a com- 
plex quantity). For such a quantum oscillator, where the 
variable c q (f) plays the role of the position variable, the 
contribution to 0S(t, q) is 



where 



pS A (t,q) = jrfiu) cos(wi), 



/M = -ft coth— -. 



(22) 



(23) 



Here, for q — > 0, lu ~ eg, and (f22|) behaves like g 2 . There- 
fore, the polaritons do not contribute to the retarded 
asymptotic form of S(t, R). 

On the contrary, in the non-retarded case c — > 00, = 



/•\/2, and (|2"2")) becomes 



PS A (t,q) = —f(uj s )cos(u) s t), 

Z7T 



(24) 



in agreement with 



IV. (cq) 2 <lu 2 <lu 2 + (cq) 2 

Now k\ is pure imaginary while k 2 is real. With q = 
(q y ,0), in region A 2 the general form of component E x 
must be 

E 2x = exp(iq y y)[a(t) cos(k 2 x) + b(t) sin(k 2 x)} + c.c. . 

(25) 

Using the Maxwell equations © , © gives the other non- 
zero components of the fields in region A 2 as 

1 k 2 

E 2y = exp(iq y y) — [— a(t) sm(k 2 x) + b(t) cos(fc 2 x)] + c.c. , 

Qy 

(26) 

B 2z = - exp(iq v y) — [d(t) cos(k 2 x) + b(t) sin(fc 2 x)] +c.c. . 
cq y 

(27) 

In region Ai , since k\ is pure imaginary, the fields depend 
on x like exp(— k\x). The continuity of E y and B z at 
x = determine the coefficients in function of a and 6, 
and ^ and (fH)|) give 



a(t) 



LO 2 - LU 2 



-4fff q ,(<), b(t) 



K\LU 

5-47rcr qw (t). 

k 2 L0 2 



(28) 

Let A 2 be the region > x > — L 2 (with a large L 2 , 
which at the end will be taken as infinite). Since U\ has 
an exponential factor exp(— 2k±x), the region Ai does 
not contribute to the total energy in this limit. The total 
energy only is F qw = A f° L dx U 2 . Using (|25H28p in $2 
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gives 



region Ai is 



iJ qw = 2irAL 2 



UJ 2 {u 2 + ~UJ 2 ){uj 2 ~U 2 _) 

uj 2 (cq) 2 \u 2 - (cq) 2 } 



kq"(*)| S 



AL 2 Cqu 



U 



2l*q<-(*)| S 



<M*)I S 



,2 \ i u 



(29) 



where w 2 is given by (JTHJ) and is given by a modified 
(fT8| with a + sign in front of the square root. Therefore, 
the contribution of this mode to (3S(t, q) is 



2^qw 



-/(w) cos(uit). 



(30) 



There are an infinity of modes of this kind, labeled by 
their frequency u, or, equivalently, by k 2 . For obtain- 
ing the total contribution of these modes, /35s (t,q) to 
/35(t, q), one has to take the sum on k 2 of (f30|) . Since 
L 2 is large, ^ fe • • • = (L 2 /n) J dk 2 ■ ■ ■ ■ Furthermore, 
dk 2 = d(uj 2 )/(2c 2 k 2 ). Therefore, 



f 



pS B (t,q) = - I d(u 2 ) 

n J(cq) 2 

x/(w) cos(iot). 



1 



n J(c q ^ ' 2cy / u 2 - (cq) 2 C qw 



(31) 



After the change of variable uj 2 — (c<7) 2 (l + u), (f3"Tj) be- 
comes, in the limit q — > 0, 

/*OC 7 

/55 B (t,q)^-^/(0)cosO / 7 =— =f-. (32) 
47T 7 s/u(l+u) Ait 

will be found the only contribution of order q to 
/35(i,q). 



V. lu 2 > uj 2 p + (cq) 2 

Now, k\ and &2 are both real. The general form of the 
components x of the electrical fields are 

Ej x — exp(iq y y)[a,j (t) cos(fcjiz:) + bj(t) sin(fcjx)] + c.c. . 

(33) 

(I3"3")) involves 4 coefficients, instead of 3 in Ei X of Sec- 
tion IV (where E\ Xl not written explicitly, involves only 
1 coefficient). Therefore the equations which were used 
in Section IV are not enough for determining all the co- 
efficients in the fields as functions of cr qw (t) . Fortunately, 
parts of the fields are uncoupled to a and, for describing 
the modes coupled to a, it is enough to choose in (|53^) 
bi = 0. Then, the Maxwell equations and the conditions 
at x — determine all the fields as functions of <r qa ,(t). 
Assuming region Aj to be of large length Li in x, one finds 
the total energy H = 7Ji qw + H 2 ^ u where the energy in 



H 



ALi W 



■<4(cqy 

= ALlClq U 

and in region A2 



UJ (uj -oj ) 



kq*(t)| 2 + -2,*q*(*), S 



(34) 



2qcj 



AL 2 C 2 qui 



|a qw (t)| 2 + -5 |<7 qa ,(t)| : 



K«(t)| 2 + ^2l*qa,(*)l' 



(35) 



For obtaining the total contribution of these modes, of 
different frequencies a;, to (3S(t, q), it would be necessary 
to sum on u the quantity 



/35 Cw (t, q) 



1 



lqu> + L 2 C 2<iuJ 



-f(u)cos(ut). (36) 



We were not able to perform this sum. Fortunately, for 
w close to Up, ([55)1 is negligible compared to (|3~4"|) and it 
will turn out that the sum on u> involves only such values 
of u in the limit of small q. Thus, we can neglect the 
term L 2 C 2aiUJ in (|3"6"|) writing 



/9Sbw(*,q) 



1 



LiC; 



-/Mcos(wt), 



(37) 



and perform the sum on uj by replacing it by an integral 
like in Section IV. Therefore the contribution of these 
modes is 



5c(i,q) 



1 



d(uj 2 ) 

W p/M cos(wi) 



A-kc* uj 2 — uj 2 — (cq) 2 uj A (uj 2 — uj 2 ) 



,(38) 



After the change of variable uj 2 = uj 2 + (cq) 2 (l + u) 
becomes, in the limit q — > 0, 



Scfoq) ~ -^f{uj p )cos(uj p t) 



du 



y/u(l+u) 



(39) 



VI. LONGITUDINAL MODES 

In addition to the transverse modes studied up to now, 
in region Ai, there are longitudinal modes of frequency 
Up (then e(u p ) = 0) @, @]. They are solutions of the 
equations (g][5|) with now p ^ and B = 0; therefore 
they occur also in the non-retarded case already studied 
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in These longitudinal modes are a superposition of 
waves with the electric field parallel to the wavenumber 
vector. Here, we summarize the calculations of 4], using 
a slightly different notation. 

Since all modes have the same frequency but differ by 
the component x of the wave vector, in ([3]) <7 qu (i) must be 
replaced by <7 q fc(t). Since B = 0, © reduces to V x E = 
0, and E is derivable from a potential. This potential is 
of the form 



the plasma — q/(47r). Things are different for the present 
retarded result, which comes from the fields in vacuum, 
with a cancellation of the contributions of the transverse 
waves and the longitudinal waves in the plasma. 

An open qualitative problem is to understand by a 
physical argument why the retarded quantum result at 
an arbitrary t is so simple, not involving H nor t. 
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<f> = [a(t) exp(ig a y) + ex.] sin(fca;); 



(40) 



this form ensures that the x component of the electric 
field 

E\ x = —k [a{t) exp(iq y y) + c.c] cos(fca;) (41) 

is non-zero at x = 0, thus is coupled by (|15j) to cr q fc(t) 
(in region A2, there are no longitudinal waves, thus E2 = 
0). (fTo]) gives —ka(t) — 47rcr q fc(t). After E 2y has been 
expressed from (|40p . one finds for the energy 



k 2 + q 2 



Kk(t)\ 2 + -o|<r qfc (t)| 5 



(42) 
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APPENDIX 

The energy density of a one-component plasma (TTJ 
can be obtained from the Maxwell equations (J3H7J) and 
the force equation (JSJ). Energy conservation can be ex- 
pressed as |5| 



After an integration on k, one finds for the contribution 
of the longitudinal modes 



J_d_ 

8tt dt 



(E 2 +B 2 ) + J • E = -V • S, 



(A.I) 



1 f K k 2 

/3S D (t, q) = ^fi^p) cos(ay) J dk p-^- 



/2 ' 



(43) 



where A" is a cut-off beyond which the use of collective 
variables breaks down. For small q, the integral in (143|) 
is K — (n/2)q + 0(q 2 ); therefore, (3SD(t,q) has a kink 
singularity at q = of the form 



where 



S = _ (E x B) 

47T 



(A.2) 



is the Poynting vector (energy flow) and J • E is the rate 
per unit volume of doing work by the electric field. ([8]) 

gives 



PS D (t,q) —f(uj p )cos(uj p t). 



(44) 



J E = 



2ird{J 2 
to 2 dt 



(A.3) 



VII. CONCLUSION 

([39]) and (|44|) cancel each other. The only contribution 
of order q to f3S(t, q) is Thus 



4tt' 



(45) 



the classical non- retarded result at t = 0, as announced 
in the Introduction. This classical non-retarded result at 
t = was found previously as the result of two contribu- 
tions: The one from the non-retarded surface plasmons 
2q/(47r) and the one from the longitudinal bulk modes in 



The lhs of (jA.ip is the time derivative of the energy den- 
sity, thus which is (fl4|) . 

The last term of (fT4| can be expressed in terms of E. 
One has to be careful about the sign: since oscillatory 
terms are to be discarded, J 2 — |J| 2 and E 2 = |E| 2 . 
Using again ([8]), at frequency u>, one can replace the last 
term in (TTJ} by ui 2 E 2 /(8itlo 2 ). Taking into account |T 
we obtain 



(ko 



(A.4) 



as written in 
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